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A Bethe Ansatz study of a self dual Zn spin model is undertaken for even spin 
system. One has to solve a coupled system of Bethe Ansatz Equations (BAE) 
involving zeroes of two families of transfer matrices. A numerical study on finite size 
lattices is done for identification of elementary excitations over the Ferromagnetic 
and Antiferromagnetic ground states. The free energies for both Ferromagnetic and 
Antiferromagnetic ground states and dispersion relation for elementary excitations 
are found. 



I. Introduction 

The present model was first proposed in 1982 by V.A. Fateev and A.B. Zamolodchikov 
[1] as a two dimensional self dual Zn lattice model with nearest neighbor spin-spin inter- 
action. Baxter, Bazhanov and Perk [2] discovered a set of functional equations involving 
families of Chiral Potts (CP) transfer matrices. Fateev Zamolodchikov model (FZM) was 
shown to be a non-chiral self dual limit of Chiral Potts [3, 4]. The Chiral Potts transfer 
matrix functional equations were used to obtain transcendental equations (Bethe Ansatz 
Equations) [5] for the zeroes of the transfer matrices for the present problem [3, 4]. In 
the general FZM model the Bethe Ansatz Equations (BAE) are coupled involving two 
automorphically connected families of transfer matrices as in the CP case [2] . In the odd 
spin case these families are connected by simple transformations and the BAEs greatly 
simplify and decouple requiring us to solve only one set of equations [4] . For this odd spin 
case alone, Albertini obtained the ferromagnetic ground state [6]. A unified treatment 
for the ground state of odd and even spin FZM can be found in a previous work [4]. 
For N = 4 FZM a comprehensive study has been done. This includes determination of 
exact energy values and central charge [7] , and completeness and classification of Bethe 
states [8]. The present work demonstrates the study of these Bethe Equations in both 
ferromagnetic and antiferromagnetic cases for finding the free energy and elementary 
excitations. 



In a generic situation for Fateev Zamolodchikov model (to be specific we shall use the 
even spin case) we obtained coupled Bethe Ansatz equations, 
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where {Xj} and {Aj} are the spectral variables for transfer matrices T q and Tjj respec- 
tively and 7 = 2lv an< ^ s ~ \ - 
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The first equation is quite similar to the generic case of Bethe Equation. However 
it involves zeroes of two different transfer matrices (A, A) that are coupled further by a 
second equation whose form is quite unique. 

The generic case of most commonly encountered Bethe Ansatz Equation looks like 

smn(Aj — ZD7J 



n sinh(A 3 - Afc - ij) = ,_, M+1 
; sinlHA, A/ • /-: 1 ' 



sinh(Aj + iSj) 



where S is the spin and 7 is the anisotropy parameter of the model. In Bethe's original 
paper [5], he studied the case 7 — > 00, where the hyperbolic functions reduce to rational 
ones. Note that in the generic case we have only one type of Xj as opposed to two species 
of Xj and Xj as they appear in the present problem. 

The standard procedure of calculating the physical quantities, e.g., the energy spec- 
trum, dispersion curves, free energy, is to assume that the solutions for the Bethe Equa- 
tion are given by the string hypothesis. Starting with the work of Bethe there has been 
a great deal of study in these complex solutions of BAE. They appear in the form, 

& = AS 1 " + \ (n + 1 - 2fc)z + + S^f k = 1, 2 . . . n 

where A^™' 1 ^ is the real part, n is its length, k runs from through n labelling the root. 
The coefficient v takes on the value (+1) (positive parity) or (—1) (negative parity); 
^oTfc^ van i snes regularly as M — > 00. 

However Bethe himself realized that for large M not all solutions of BAE are of the 
form (above) with UmM—>ooIm(5) = 0. Modern work for the case of limM—>ooIm{$) 
was initiated by Destri and Lowenstein and Woynarovich who introduced the definition 
of narrow pairs for lim m-> 00 1 m(S) < and wide pairs for limM^ooIm(S) > and was 
furthered by Avdeev and Dorfel [9, 10] who introduced 3 classes for lim m-> 00 1 'm(§). It 
is clear that none of the existing approaches to BAE for S integer or half integer is suf- 
ficiently refined to answer the reality or completeness question for the Hamiltonian. 

In the general spin hypothesis framework one obtains equations for the centers of the 
strings by multiplying the Bethe Ansatz equations over different members of the same 
string and then taking the logarithm of the resulting equation. This yields 

1 1 M<fc) T U) 

k 13=1 

where is the number of fc-strings and 

rij 

&f ] (X) = 2J2 <t>{\ rij+2s-2l+l, V j) 

minOj ,n k ) — l 

©jfc(A) = <t>{X,rij +nk,VjVk) + <A(A, \rij - n k \,VjV k ) + 2<f>(\, \rij - n k \ + 2l,VjV k ) 

and 

4>(\,n,v) = l 2^rctan(cot(^)nanh(A)) 
v ' if ny = qn q e Z 
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In certain cases, e.g., <5-function Bose gas, it can be proved that the solutions are real, and 
no such multiplying of string components is necessary [11]. In such cases, the integers des- 
ignating the branches of logarithms uniquely characterize the states and may be viewed 
as quantum numbers for the states. A monotonic relation is shown to exist between the 
integers and the values of the spectral variable a£™'^ . In almost all subsequent work in 
the field this unique characterization of states by integers and their monotonic relation 
to solutions have been assumed. In few cases some counting argument is attempted 
[5, 12, 13, 14, 15] to justify this assumption. In the N — > oo limit, after introducing 
the concept of density of string centers, one obtains a coupled set of integral equations. 
These equations are manipulated to calculate the energies of the ground state and low 
lying excited states. For the spin 4 system, complete classification of Bethe states, and 
exact calculation for finite and infinite systems is already known [7, 8]. 

In the present problem, we had to deal with a doubly coupled set of integral equa- 
tions; two sets of coupled equations both involving zeroes of two types of transfer ma- 
trices. However linearity of the equations made it possible to solve for the ground state 
and elementary excitations by Fourier transform method. A study of the classification 
of roots is undertaken. However one has to keep track of the added complexity of having 
to handle T q and Tjj simultaneously. One can identify the ground state and elementary 
excitations on the basis of this numerical study. The excitation spectrum and dispersion 
relations can hence be calculated. 



II. Fateev-Zamolodchikov model 

V.A. Fateev and A.B. Zamolodchikov proposed in 1982 a two dimensional self dual Zm 
lattice spin model with nearest neighbor interaction. They obtained this model as the 
self dual [16] solution of the star-triangle relations [17]. 

A general model can be defined as follows. On a two dimensional rectangular 
lattice the lattice sites are occupied by a spin variable z which takes its values in the 
group Zn [z n = 1]. If one designates the sites on the lattice by a two dimensional 
integer-valued vector x, one can write down the partition function of the statistical Zn 
model with nearest neighbor interaction as: 

z =Enn ^w^)'^^))- ( 3 ) 

{z} X <j=± 

where the sum runs over all values of the variable z in every site of the lattice. The 
functions w a , (a = ±1) are the weight functions corresponding to the interaction 
between spins on the neighboring sites of the lattice in horizontal (a = 1) and vertical 
(a = —1) directions respectively. The vectors ei = (1,0) and e_i = (0,1) are the basis 
vectors of the lattice. 



In the absence of external fields, the most general interaction between two neighboring 
spins after appropriate normalization is given by 

N-l 

w^{z 1 ,z 2 ) = l+Y^x { ° ) -{z.ziy (4) 

where superscript * denotes complex conjugate. Reality of w^(zi,z 2 ) imposes on the 
parameters the following restriction 
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Figure 1: Fateev Zamolodchikov Model on a square lattice 
The dual transformation of the statistical weights are given by 



~(<7) 



N-l \ / JV-1 \ 

fe=l / V k=l / 



where w = exp(2iri / N) . The region of self duality is then given by 



x^=x^ 



(6) 



(7) 



Let the parameters x^ be represented by a family of functions Wi (a) of auxiliary pa- 
rameter a G C 

x [ p = Wi(a), x { r^ = W t (ir - a). (8) 
The star-triangle relation [18, 19] on Xj(a) : 



N-l 



J2 W ni - k (a)W n2 - k (TT -a- a')W n3 . k (a') 



k=0 



c(a,a')W n2 - n3 (n - a)W ni - n3 (a + a')W ni - n2 (ir - a') 



(9) 



The particular solution of Eq. (9) that possesses the self duality property e.g. Eq. (7) 
is given by: 

sin[7rfc/iV + a/27V] 



W = l, W n (a) = T[- 



sin[7r(fc + l)/Ar-a/2AT]' 



Denoting x„ ^ — W(n \ u) and x [ n L> = W{n \ u) we get 

W(n\u) _ -A- sin(7T j/N - ir/2N - it) 
W(0\u) ~ AA sm(nj/N - 7r/2iV + it) 



(-i) _■ 



W{n\u) _ -A- sin(7rj/7V - n/N + u) 

w(o\u) A! 



sm(n j/N — u) 



(10) 

(11) 
(12) 



We adopt the normalization Vt^(0|u) = VF(0|it) = 1. The "physical region" defined by 
non-negative real Boltzmann Weights (BW) , is given by u e [0, tt/2N[. For N = 2, 3 Eq. 
(11) and Eq. (12) simply reduce to the self-dual critical Potts model. For N = 4 it gives 
a particular case of critical Ashkin- Teller model. Fateev and Zamolodchikov propose that 
for N = 5, 7 the solution describes the critical bifurcation points in the phase diagram 
of Alcaraz and Koberlc [20]. 
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Figure 2: Square lattice Chiral Potts Model 



III. Chiral Potts model and connection to FZM 

On the sites of a two dimensional lattice of size M x j\f denoted by two dimensional 
vector (j, k) with integer entries, we place Zjy spins aj t k- The spins aj^ are classical 
variables satisfying: 

a 3-k = 1 

i.e., Cj.fc = lo v ', ^ G {0, 1, 2, ... JV — 1} where w is the complex N th root of unity with the 
minimum argument. 

The energy corresponding to a given configuration of spins {<Jj,k} is 

JV-l 

£ = ~ E E ^ ■ K^^+l)" + ^ • (^fc^l,*)"}- 
{j,fc} n=l 

Row index j runs over 1 to M and column index k runs over 1 to TV with periodic 
boundary condition in both directions implied. 

In Chiral Potts a subspace of the coupling parameters (E%, E£) is chosen which has 
a built-in handedness, or phase. The energy of a nearest neighbor pair is chosen as 

n=l 

where E% v = \E*' V \- e i<5 ". The local Boltzmann weights (BW) can now be easily defined 
as: 

W h ' v (n) = e*5r^J=i ' 

Let us denote two adjacent row configurations by {1} and {/'} ( I corresponds to the 
lower row), where 

{1} = {Ji | j = 1(1)tV and lj G 0, 1, . . . , N - 1}. 

The row to row transfer matrix is given by: 

TV 

T{i},m = ]lW h (l j -l' j+1 )-W v (l j -fy. 
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Figure 3: Transfer Matrix 



It has been shown by several authors [1, 21, 22, 23] that the transfer matrices corre- 
sponding to the interaction parameters belonging to the Chiral Potts submanifold com- 
mute. 

[T, T'} = 

The self dual Chiral Potts model is given by BWs 

Wpgjn) = tt b q -^a p 

W pq (0) W&p-u^V 1 j 

^ggW = TT ua p -u3g q 

w pq (0) & 9 - ' 

where o> = cxp(27ri/7V) and the paired complex variables (a, 6) £ C 2 satisfy the con- 
straint 

a? + b» = K (15) 

k G [0, 1], and x = p or q. In the non chiral limit, when n = 0, we can parametrize 
(a x ,6 x ) in Eq. (15) as: 

a x =e 2lx b x = uj 1 ' 2 e 2lx . (16) 

Defining u = q— p Eq. (13) and Eq. (14) reduce to Eq. (11) and Eq. (12). However 
we retain suffixes (p, q) in the Boltzmann weights W pq (n\u) and W pq (n\u) to signify that 
these BWs are obtained from the Chiral Potts BWs defined in terms of p and q variables. 
The transfer matrix for the FZM can be constructed from the BWs as: 

M 

T ^ n \ u ) = l[W pq (n k -n' k \u)W pq (n k -n' k+1 \u). (17) 
fe=i 

where M is the number of sites in each row and periodic boundary condition is implied. 
These transfer matrices for different spectral variable u form a commuting set. This can 
be argued from the fact that these transfer matrices come as a limit of CP transfer matrix, 
which are known to be commuting. A more direct argument would be that Fateev and 
Zamolodchikov obtained FZM BWs as solutions of star-triangle relation, and hence the 
transfer matrix constructed out of them ought to commute. 

[T(u),T(u')] =0 V u,u' G C (18) 



G 



Transfer matrix T(u) reduces to identity operator for u — > 0. An expansion of T(u) 
gives us the associated spin chain Hamiltonian H. 



N-l 



T{u) = 1 - Mu Y . , 1 — - uH + 0{u 2 ) (19) 
^ sin(n7r/iV) 

n=l K 1 ' 



M N-l 



fe=l n=l v ' 7 



where X and Z are defined as 



X k \m . . . n k . . . n M > — \n\ . . . n k + 1 . . . n M > modiV 
Z k \m . . . n k . . . n M > = uJ nk \n 1 . . .n k . . .n M > ■ 



Eq. (19) and Eq. (20) imply that each Hamiltonian commutes with all the transfer 
matrices and their associated Hamiltonians. Thus it has an infinite set of conserved 
charges in involution. However only a subset of them, whose number is equal to the 
degrees of freedom of the system, are independent. 

In order to obtain the zeros of the eigenvalues of the transfer matrix T q , we will use 
functional equations connecting T q with its automorphically conjugate partners. Thus 
it is important to understand the relevant automorphisms of the constraint Eq. (15). 
It has been claimed in the previous section that the transfer matrices constructed out 
of CP BWs, Eq. (13) and Eq. (14). commute as long as they satisfy Eq. (15). For 
any (a, b) <G C 2 satisfying the above relations there exist other complex pairs connected 
to them which satisfy the same relation. Two such automorphic relations of importance 
are, 

R{a, b) = (6, ua), U(a, b) = (uja, b) (21) 

It is rather straightforward to check; 

n , h N JV I — 

a Rx + °Rx — K i a Ux + °Ux — K 

from the relation 

+ b% = k 

If one makes an attempt to go over from CP BWs to FZM BWs through a limiting 
process, one gets the following relations for W P R q , W P R q , W p u q and W p u q , 

W pRq (n\u) _ -A- sin(7rfc/7V - ir/N - u) 



^(Trh IN 4- n\ ^ ' 



W pRq (0\u) 1JL sm(7rk/N + u) 



Wp Rq (n\u) yr sin(7rfc/jV - 7T/2N + u) 

W pRq (0\u) ~ l = lsin(7rfe/JV-7r/2JV-u) 1 ' 

W pUq (n\u) ^ " S m(irk/N-Tr/2N-u) 

W pUq {0\u) ' ' l} iS m(TTk/N + ir/2N + u) [ ' 

Wpu q (n\u) eW -A- sin(7rfc/jV + u) 

W pUq (0\u) " sin(7rfc/7V - u) 

Thus in the non chiral limit, T q — ► T q {u) and T Rq — > T q (u + n/2N). There is no simple 
relation between T q and Tjj q though. However, we do feel that there must exist some 
nontrivial mapping between T q and Tjj q whose understanding will unravel the connection 
between the zeroes of T q and Tjj q and will give the satisfactory derivation of completeness 
of states. 
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IV. Bethe Ansatz type equations for the even N FZM 

We define the normalized transfer matrices by removing their denominators, 



T q N (u) = [g q (u)g q (u)] M T q (u) 



where 



N/2 . N/2 

3=1 3=1 



One must note that the superscript in T q denotes "normalize" and is not related to the 
spin quantum number N. Each entry of T q (u) is a product of NM sines and it has the 
general form 



NM 

n< 

k=i 



The calculation of this section goes in the same spirit as that of odd N case. Hence 
we only quote the results. 



Aq=o(u) = 



,(0)5,(0) 



M L 



g q (u)g q (u) 

L = A+B= NM 

The normalization has been fixed by T q (0) = 1^. 
The momentum (P) is given by, 



n 

k=i 



sin Vk 



(26) 



iP 



2N> 



-.ML 



9q( 27v)f<?(27v) 



n 

fc=i 



^2N 

sin Vk 



(27) 



Now we turn to the sectors Q ^ 0, and the symmetry under charge conjugation allows 
us to consider the sectors Q = 1,2, . . . ,(N — l)/2 only. While we have not been able to 
obtain a proof like the one given above, one can show that, in the sector Q 



a) A, B < 

b) A, B > 



NM 
2 

NM 



Q 



^■4 



N 
~2 



Following similar argument as before, we arrive at 

NM 



A = B = 



-Q 



The reader must be warned that this conclusion lacks rigor just like in the case of N 
odd. The factorization in terms of sines can be carried out without the appearance of a 
phase ( e 2 »«)±(B--4)_ ^y e assume this to be true also for the others Q sectors, and arrive 
at the general form 



Aq(«) 



,(0)9,(0) 



M L 



n 

fc=i 



sin (it — Vk) 



_g q (u)g q {u) 
L = NM -2Q, Q = 0,1,.. 



sin Vk 
N 



(28) 



Ajv-q(w) = A Q (u) 
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From this, the eigenvalue of H is easily found to be 



L N/2 

E = ^cotv k - 2M'Y^cot{iTj/N) 

fe=l j=l 



(29) 



The momentum (P) is given by, 



2N > 



M L 



n 

fe=i 



< 2N 

sin Wfe 



(30) 



We shall use the set of functional equations for the eigenvalues of transfer matrices 
of Chiral Potts derived by Baxter, Bazhanov and Perk [2, 24]. This functional relation 
appears in reference [2] as Eq. (4.40) and has the following form, 



N-l 



T q — c m ^ q T umq T q T Um+lq X 



m — l 



(31) 



m=0 



where f — TS, q — (a,, bq) ~UR 1 (a q , b q ), and 



M 



j=0 



, j=m+l 



b v — uji +l a a 



a p b p — uj m a q b q 



2 2 



s=0 s=0 

2 1 / 

1/4 A , R 23 q 1 1 R 2 °{Uq)^ + A ,A ' 1 R?°(Uq) 1 1 1 R 2 °+ 2 q A 

s=0 X^-s.qSis.q' S±s,q' S±s+l,q 



Define 



P2s 



C2s,q 



A A , 

^-s^q^s^q' 



d: 



2s+l 



c 2s+l,q 



A s ,q' A s +l t q 

The independent inverse factors of T R 2, q and T R 2B^ Uq ^ are considered, and both sides 
of the above equation are multiplied by the appropriate common factor so as to get rid 
of inverses of transfer matrix. The appropriate factor is: 



(32) 



After multiplying we get 



X m Tq ■ T R 2 jq T R 2(j-i)( Uq ) — ^2 
j=l s=0 



v 



3 = 



3 = 



\ 



Cfes + l { [ T R 2 Jq ■ Y[ T R 2j 



(Uq) 



3 = 

j¥(*+i) 



3 = 



(33) 
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If one expresses T v>q and T v jj q in terms of a complex parameter u, where u = q — p as 
fc^O, 



T q -> T q T Uq -> T Uq (u) 

TV 



T fl 2 N T 9 (m + ^) T R2 k {Uq) -> T (;7g) (u + ^) (34) 



with this parametrization we get 



1 f-i / 4 4- 



X.f q .]jT q (u+^).l[T Uq (u+^)=Y^ 

j=l j=0 s=0 



^n^+f ) n 

3=0 j = 



V 



f . f-1 

+rf 2s+1 n t,( U +^) n i (35) 

3 =0 3=0 

Let w be a zero of T q , i.e., T 9 (w) = 0, then whenever u = v — irk/N, T q (u + nk/N) = 0. 
Thus for u = v — nk/N k G {1, 2, . . . , ^ — 1} all but two terms vanish. 

JV JV_ JV N__ 

3 = 3=0 3 = 3 = 

Cancelling the common factors we get 

7T k TT IT k 

P2k{u) ■ T Uq (u + — - -) + dafc-^u) ■ 7V,(u + — ) = (37) 



whence 

T Uq (v) P2k(v-f) 



T Ug (v-§) d 2k ^(v-f) 
Recalling the expression for Tu q , 

sm(vi-Vj) _ ( 9 Uq ( v i) ' 9uA v i) \ M P2k(v l -^) 



(38) 



n 



fJi sin («i - vj - f ) \g„M - #) • g Uq (vi -jf)J d 2k -i(vi - 
The ratio of g Uq 's can be obtained as 



(39) 



V i V o 

2 - 2 - 

/ \ TT . /KJ ft \ / TT \ TT • fftJ 71 \ 



N 2N in * Uqy 1 N J- 1 ./V 2AT 

_ sin (§ - f + m + v i) _ cos ( w * - 



ff^K-f) sin(0 + 2^ + Ui) sin(wi+2^) 

Similarly the ratios of g Uq is found as, 

N_ -1 JV 

2 1 2 

a uq m = n sin (77 - 9 Uq ^ - = n sm (77 - «o 

ff„ 4 («i) _sin(f-u i ) sin(« i + $) 



(40) 



^(^i-f) sin(§-^) cos(«i) 



(41) 
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Using these results for the ratios of gu q and those of (p2k/d2k-i) we finally get the Bethe 
equations [5] 



n 



sin(ui - Vj) 



(-1) 



M+l 



- sin2fa-^) 
sin(2wi) 



2M 



(42) 



Let v be a zero of T Uq , i.e., T Uq (v) = 0, then whenever u = v — nk/N, T q (u + nk/N) = 0. 
For u — v — nk/N k € {1, 2, . . . , ^ — 1} all but two terms vanish 

JV — — 1 — — — 1 

P2 k f[T q (u+^) [] T c/9 ( U +^) + d 2fe+1 n T q (u+^) n 7^+^ = (43) 



TV' -LJ- TV' 

3 = 3 = 

Cancelling the common factors we get 

irk 7T 



3 = 
jjik+1 



3 = 



P2k{u) ■ T q (u + — + -) + d 2k+1 (u) ■ T q (u + — ) = 
T g (v) 



T q (v + f) " d 2 k+1 (v-f) 



(44) 



The ratios of g q 's can be obtained as 



7,(150 = 11 sm(— - — + Vi), 



3 = 1 



9 v (vi) sin( 



jV 2N 



9 q (v l + -)= 11 sm(- - — +«,) 

j=2 



+ «») sin(^+2^) 



g q (vi + §) sin(f + f + t;,) cosfa + ^) 
Similarly the ratios of g 's is found as, 



(45) 



sin(f-Ui) 



9„( v i + 7i) sin(0-w,) 



2 

- /- ^ \ TT • /^J _ 



(_1) . — ^ 



N 



sia(vi) 



(46) 



Using these results for the ratios of g q and those of P2fc/^2fc+i we finally get 

yr Sin(f; t - Vj) = M+1 

11 sin( ^. _ v . + * ) ^ ) 

In order to cast the BAE's for even case in a simpler (and standard) form, we make 
a change of variables. 



The BAE's in terms of these new variables are 



Vj — i\j — — 
3 3 4N 



n 



sinh(Aj — Afe — 27) 
sinh(A j - A fe + ij) 

yr sinh(A 3 - Afe - 27) 
sinh(A j - Afe + i 7 ) 



(-1) M+1 
(-1) M+1 



sinh2(Aj + is7) 
sinh2(Aj — is7) 



2il/ 



(47) 

(48) 
(49) 
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where 7=2^ an d s = \- From the numerical study for even spin BAE-s, it was found 
that XjS are related to one another. In fact 

V Xj 3 Xj+iir/2 mod(Tr) (50) 

This allows us to group Xj such that Xj G [— 7r/4, 7r/4]. Using transformation rules for the 
hyperbolic functions one can rewrite the expressions in terms of a new variable Xj = 2Aj. 
The LHS of BAE(l) becomes 



sinh(Aj — Afe — 17) sinh(Aj — Afe — 17 — ^ ' 



n 



^ sinh(Aj - A fe + 17) sinh(Aj - A fe + ij - *f ) 



= -r-r sinh2(A 3 - - A fc - ij) = -pr sinh(xj - Xfc - gry) 
11 sinh 2{Xj - X k + * 7 ) W sinh( Xj - X fc + 2* 7 ) 

RHS of BAE(l) is rewritten in terms of variables Xj 

M+i ( sinh(x 3 +2is7) \ 2M _ M+1 / sinh( Xj -+2 7 ) \ 2M _ 1 

A similar transformation is done for BAE(2). Hence the BAE equations become 

ff sinMxj - Xk - 2^7) = , ^M+i / sinh(xj +H) \ 2M ^ 
W sinh( Xj - Xfc + 2i 7 ) ~ ^sinh( Xj -i7)y 

n sinhfe - X fe - gyy) , ,m+i 

,sinh;\, u • ■>,-; [ > [ > 



The eigenvalues Aq(u) of the transfer matrix are given in terms of these zeroes Xj s 
as [3], 



Aq(«) 



g g (0)g,(0) 

9q{u)9 q {u) 



M 



n 



sin(2u - ^Xfc ~ 2^f) 
sin(iXfc + w) 



(55) 



where, 



N/2 

9q( u ) = II sm ( 



7TJ 7T 

77 ~ 2A 



JV/2 



(56) 



V. Study of finite size systems 

The BAE in the present model differs from the standard form. The sign in front of 7 on 
the right hand side (RHS) of BAE (1) is reversed. That is to say that the RHS is equal 
to the inverse of what usually is known to be the RHS. The second of the coupled pair, 
BAE (2) is even more striking. Though the LHS is still coupled, the RHS is independent 
of the spectral variable ! Understanding the real significance of these peculiarities can 
help enormously in solving the problem. 
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The transfer matrix for the FZM is constructed from the FZM Boltzmann weights 
(BW) as: 

M 

If '"'(«) = T^ n '(u) = l[W pg (n k - n' k \u)W pq (n k n' k+1 \u). (57) 

k=l 

where M is the number of sites in each row and periodic boundary condition is implied. 
These transfer matrices for different spectral variable u form a commuting family. Trans- 
fer matrix acts on vectors defined in terms of spin indices along a row (or diagonal) [3] 
or the spin configuration n = |m, n 2 , ■ ■ ■ nu >■ 

There is an associated transfer matrix T p ,u q which corresponds to a conjugate set of 
Boltzmann weights [4], 

/ / M _ 

Tu q ' n («) - (u) = [] W pUq (n k - n' k \u)W pUq (n k - n' k+l \u). (58) 

fe=i 

The shift or translation operator S is defined by its action on a state function or spin 
configuration n = |m, n 2 , ■ ■ ■ Hm >■ 

S\ni,n 2 ,n 3 , ...n M >= \n M ,n 1 ,n 2 , ■ ■ .n M -i > (59) 
Momentum P is defined in terms of the shift operator as 

e iP = S- 1 

The z-componcnt of spin operator Z k and spin raising operator Xk corresponding to 
a given lattice site (fc) are defined by their action on a state function or spin configuration 
as n = |m, n 2 , . . . tim >, n k = 0, 1 . . . N — 1. 

X k \m . . . rife . . . riM > = \n\ . . .n k + I . . .hm > mod(A f ) 
Z k |m ... rife . . . n M > = u nk \ni ... n k ... n M > 

The global spin raising operator is given by 

M 

x = ]Jx k 

k=l 

The spin Q is defined in terms of X as 

= X- 1 

T p , q (u) and T Pt u q (u) commute with S, Z k and X. Hence it is possible to make a spin 
and momentum sectorwise study of the problem. 

The roots of the BAE are studied by computing the eigenvalues of the transfer ma- 
trices as meromorphic functions of x = e x . Since the transfer matrices with different 
spectral parameter A commute, their eigenvectors arc independent of the spectral pa- 
rameter. Hence by taking a specific value of the spectral parameter one can determine 
the eigenvectors numerically by diagonalizing the finite size transfer matrix. From the 
definition of the eigenvalue equation one can express the eigenvalues as meromorphic 
functions of x, as the entries of transfer matrix are polynomials in x and the eigenvectors 
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are vectors with numerical (independent of x) entries. 

However the problem being coupled, one needs to simultaneously diagonalize T q and 
Tu q or in other words find the eigenvalues corresponding to simultaneous eigenvectors of 
T q and Tjj q . This fact by itself introduces a significantly higher level of difficulty over 
other numerical simulation of similar type (non coupled eqns.) e.g. Chiral Potts [25]. 
Coupled BAE and any simultaneous eigenvalue problem (eigenvalues corresponding to 
simultaneous eigenvectors) results in the same generic problem in computer algorithm. 
One has to develop efficient optimized codes for tackling this. 

A detailed numerical study for chains of length M < 8 [3] was done. The main 
observations of this numerical study are as follows: 

• 1-String with both parities: (l,v) ,v = (±1) 

• even length strings with positive parities: (n, +), n = 2, 4..., N 

• non-string solutions Im{\) ~ ±7r/3 

The Ferromagnetic ground state is a filled band of 2-strings, and the excitations con- 
sist of (1+), (1 — ) etc. The Antifcrromagnetic ground state on the other hand is a filled 
band of (1, ±) and excitations are 2-strings with positive parity. 



« 


2 

> < 


-strings 
» (2s) 


i 


^ non-string 
(ns) 


• 

real -solution 
(1+) 

< 


t i 

(l-)s 


» 

tring ^ 


• 

« 





Figure 4: String and non-string solutions for BAE 



It is remarkable that a good deal of insight into the nature of string solutions can 
be obtained from finite systems of rather small size. The 2-strings are easily identified 
as having imaginary parts close to 7r/4. One can also identify the real roots and roots 
with negative parity Im{\) — n/2. There are roots whose imaginary parts are not well 
approximated by ir/2 or 7r/4. These roots do not seem to systematically approach the 
2-string value as M — > oo. These roots are classified as non-strings. For space limitation 
M = 2 and M = 3 figures are presented in this paper. The imaginary parts of the 
non-string (ns) roots are 

M=2 P=0 -0.32198235™ 

0.321982357ri 

M=3 P=0 -0.297309027ri 

0.297309027ri 
P=l 0.179912 - 0.3305457ri 

0.179912 + 0.33054515?ri 
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P=2 -0.179912 - 0.330545™ 

-0.179912 + 0.33054515™ 
M=4 P=0 -0.33780988™ 

0.33780988™ 
-0.348355 - 0.342016™ 
-0.348355 + 0.34201621™ 
0.348355 - 0.342016™ 
0.348355 + 0.34201621™ 
-0.28324874™ 
0.28324874™ 
P=l -0.258619 - 0.38794057™ 

-0.258619 + 0.38794057™ 

The following tables show the classification of roots for spin-4, for lattice sizes 2, 3 
and 4 in the Q = sector. The first column shows the momentum P. The second column 
gives the roots of the BAE in the variable \- The content or the type of root is identified 
in column three. Column four gives the corresponding integer that appears in the BAE. 
In column five the energy calculated for the corresponding eigenvalue for Transfer Matrix 
is given. 



Table 1: Classification of roots and integers for M = 2 



p 


A fe = lnx fe 


Content 


h 


Energy 




-0.216337- 0.227395™ 


(2s) 


0.5 






-0.216337 + 0.22739469th 










0.216337 - 0.227395™ 


(2s) 


-0.5 







0.216337+ 0.22739469™ 






-6.24622 




-0.265319 


(1+) 


-0.5 






-0.32198235™ 


(ns) 









0.32198235™ 











0.265319 


(1+) 


0.5 


4. 




-0.440687 


(1+) 


-1. 









(1+) 









— 0.57ri 









1 


0.440687 


(1+) 


1. 


8. 



Table 2: Classification of roots and integers for M = 3 



p 


Afe = lnxfe 


Content 


h 


Energy 




-0.325636 - 0.2288977ri 


(2s) 


1. 






-0.325636 + 0.228896647ri 










-0.239862447ri 


(2s) 









0.239862447ri 










0.325636 - 0.2288977ri 


(2s) 


-1. 







0.325636 + 0.228896647ri 






-8.36945 




-0.388991 


(1+) 


-1. 






-0.2677187ri 


(2s) 









0.2677187ri 










-0.297309027ri 


(ns) 









0.297309027ri 











0.388991 


(1+) 


1. 


-1.65686 
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Table 2 Contd.: Classification of roots and integers for M = 3 



p 


A fe = lnx fe 


Content 


h 


Energy 




-0.570869 


(1+) 


-1.5 






-0.213538- 0.5™ 











-0.149008 


(1+) 


-0.5 






0.0873277 


(1+) 


0.5 






0.423044 - 0.242196^ 


(2s) 


-1.5 







0.423044 + 0.24219619™ 






8. 




-0.768026 


(1+) 


-2. 






-0.241758 


(1+) 


-1. 






-0.24949865™ 


(2s) 









0.24949865™ 










0.241758 


(1+) 


1. 







0.768026 


(1+) 


2. 


4.48683 




-0.353241 - 0.236642™ 


(2s) 


1. 






-0.353241 + 0.23664231™ 










-0.0132479 


(1+) 









0.179912 - 0.330545™ 


(ns) 









0.179912 + 0.33054515™ 








1 


0.359907 


(1+) 


1. 


3.37155 




-0.540092 


(1+) 


-1.5 






-0.129032 


(1+) 


-0.5 






-0.0725776- 0.5™ 











0.117714-0.250431™ 


(2s) 


-0.5 






0.117714 + 0.25043057™ 








1 


0.506275 


(1+) 


1.5 


3.10102 




-0.692847 


(1+) 


-2. 






-0.209487- 0.250419™ 


(2s) 


1. 






-0.209487 + 0.25041875™ 










0.0217873 


(1+) 









0.2683 


(1+) 


1. 




1 


0.821734 


(1+) 


2. 


6.82843 




-0.856212 


(1+) 


-2. 






-0.287304 


(1+) 


-1. 






-0.0386368 


(1+) 









0.185678 


(1+) 


1. 






0.375937- 0.5™ 




-0.5 




1 


0.620538 


(1+) 


2. 


12.899 




-0.885356 


(1+) 


-2. 






-0.304098 


(1+) 


-1. 






-0.0535889 


(1+) 









0.16349 


(1+) 


1. 






0.539776 - 0.247284™ 


(2s) 


-1. 




1 


0.539776 + 0.24728422™ 






11.4569 




-0.359907 


(1+) 


-1. 






-0.179912- 0.330545™ 


(ns) 









-0.179912 + 0.33054515™ 










0.0132479 


(1+) 









0.353241 - 0.236642™ 


(2s) 


-1. 




2 


0.353241 + 0.23664231™ 






3.37155 
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Table 2 Contd.: Classification of roots and integers for M = 3 



p 


A fe = lnx fe 


Content 


h 


Energy 




-0.506275 


(1+) 


-1.5 






-0.117714-0.250431™ 


(2s) 


0.5 






-0.117714 + 0.25043057™ 










0.0725776 - 0.5m 











0.129032 


(1+) 


0.5 




2 


0.540092 


(1+) 


1.5 


3.10102 




-0.539 776 - 0.247284™ 


(2s) 


1. 






-0.539 776 + 0.24728422™ 










-0.16349 


(1+) 


-1. 






0.0535889 


(1+) 









0.304098 


(1+) 


1. 




2 


0.885356 


(1+) 


2. 


11.4569 




-0.620538 


(1+) 


-2. 






-0.375937- 0.5™ 




0.5 






-0.185678 


(1+) 


-1. 






0.0386368 


(1+) 









0.287304 


(1+) 


1. 




2 


0.856212 


(1+) 


2. 


12.899 




-0.821734 


(1+) 


-2. 






-0.2683 


(1+) 


-1. 






-0.0217873 


(1+) 









0.209487 - 0.250419™ 


(2s) 


-1. 






0.209487 + 0.25041875™ 








2 


0.692847 


(1+) 


2. 


6.82843 



VI. Free Energy in the Ferromagnetic case for N even 

From the numerical study one can identify that the Ferromagnetic (FM) ground state 
corresponds to a filled band of N/2 strings of positive parity for T q and a filled band of 
1-string of negative parity for T\j q . This vector always falls in the P = sector as is 
expected. A further study upto 6 sites reveals that this remains true. 



The LHS of the first of Bethe Ansatz Equations, BAE(l), is given by: 



sinh(x 3 -Xk~ gry) 
LA sinh(xj - Xk + 2ry) 



n 



M+l 



sinh(xj + ji) 



-, 2M 



smh(xj - 7*) 



(60) 



For the Ferromagnetic case we made the assumption that the ground state corresponds 
to A^/2 strings with positive parity for T q and 1-strings of negative parity for Tjj q . 



Aa,v 

Y™. 1 
Aa.v 



xS + 27(n+l-2i)i 



y 

x a 2 



(61) 
(62) 



Define x and x by x — an d X = 2jx. If denotes the number of n-strings, we 

get for the left hand side (LHS) of the BAE(l) 



n 



sinh 27(2:™ — xp + (n + 1 — 2l)i — i — poi) sinh 2j(x™ — xp — ni — poi) 



^ sinh 27(^2 — xp + (n + 1 — 2l)i + i — poi) sinh 27(0;™ — xp + ni — p a i) 
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Taking product over the string elements of \a the first BAE becomes decoupled and 
is given in terms of the variables for T q alone. 



(_l)(M+l)n 



n 



sinh2 7 (.T« + (n+l-2Z + i)z) 



2M 



= (-1) 



M+l 



sinh2 7 (a;S - (n + 1 - 21 + 
After multiplying for the elements of a string, the second BAE becomes 



(63) 



nn 

*:=i i=i 



sinh 2 7 (a: a 



(n + 1 — 2Z)i — i — poi) 



sinh 2"/(x a — Xg — (n + 1 — 2Z)i + i — po*) 



(-1) 



(M+l) + - 



(64) 



Thus we have only one set of BAE which involves the zeros of the transfer matrix T q . 
From the above equation taking natural logarithm of both sides we get 



2M i ■ In 



1=1 



sinh( X a + 2 7 (n+l-2Z + i)) 



v sinh(xa - a 7 (n + 1 - 21 + |)) 

Defining the density of string centers for the zeros of T q by 

1 



P(x) = l im 



we get, 



where 



M^oo M(xfe+1 - Xk) 



p(x)=>g; +) (x) 



(65) 
(66) 



n , 

e (~, +) (x) = E 2 ^.»+2- 2i ' + ) 

and prime on 0^ ,Ax) denotes differentiation with respect to the variable \. Here the 
function 0, as defined by Takahashi and Suzuki [26] is, 

sinh 2 7 (x + to + po*) 
g(x,n,v) = -r 



<j)(x,n,v) =i-]n(g(x,n,v)) ; 
Evaluating the sum over I in the Fourier space, 



sinh 2 7 (a; — ni + poi) 



Qfl.Ak) = 2ir 



sinh(f - $) 



sinh(f) 



'(f,+) 

By inverse Fourier transform we get 

1 r +oc 

P(X) = ^ J dke-^p(k) 



N ■ 

sinh(f - f|) 
sinh(f ) 



(67) 
(68) 



The free energy for the Ferromagnetic ground state is defined as 



sin(2u — + ^(0) 



i f=i sin (*Xa - x(l)) 
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Replacing the summation by an integral over the symmetrically placed string centers, 
r, s 1 , , /cosh(2x) -cos(4u + 2x(/))\ 

Transforming to the Fourier space and utilizing the expression and properties of p(k) 
we get 



I 4 dfc sinh(fc7r - ^1) sinh(2fc?i) sinh(2fcu - kir - !f) 
J " [ " ] ~ ! ~ sinh 2 (^) ~ " "° 



VII. Excitation on FM ground state 

We have seen in the previous section that the FM ground state is given by a filled band 
of (2s) strings. Consider the Z^ 2a )(x) function. In the general case it should look like 

%*)(x) = ^eg>>(x) -^EE e$,k)<x-x&) (7i) 

k 0=1 

The density of (2s) vacancies is given by 

t { 2s)(x) = -Z' {2s) (x) (72) 
The vacancy density and the density of (2s) particles is related by 

M< 2s > 

a (2s) (X) = P(2.) (x) + ^ E «(x- xf )h ) (73) 

/3=1 

where are the position of the holes. 

Thus 

M< fc > 

-WX) = ^<)W-^M E E<, fe )(X-X^)- 

fc#(2s) /3=1 



M< 23) 



h=l 

The above equation can be interpreted as a collection of terms contributing to (2s)- 

(o) , (h) T (j) whcre(J (o) 

(2s) + a {2s) + a (2s) wnere <J ( 2s ) 



ground state, (2s)- holes and excited particles <?(2s) = a ( 2 l) + a (2s) + where cr'- ' 1 



is the same as /9( 2s ) of the last section. 

The expressions for energy (E) and momentum (P) are 



_L JV 

E = Ecot(% + ^)-2MEcot(^) (75) 

I 1 Z 1 V / 



fe=i fe=i 
sinh (Xfe + ^) 



n 2jv (76) 

l = l s inh( Xfe -^) 
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We obtain the energy of a state designated by a given set of strings 

M (k) 



0=1 

strings 



/ 



0=1 

M< 23) M ( fe ) 



^( 2s )(x)e (25 )(x) - E (4 2S)/l ) + EE e (*?>) (77) 



0=1 /3#fe /3=1 

The bare energies for n-string with parity w is easily obtained 



Hn,v) (Xa) = E COt + ^) 

fe=l 
n 

= E cot {^i n ' v) 2 7 (n + 1 - 2k) - \{l - v) + ^) 

fe=l 

n 

= E tan (f - ^"'^ + 2 ^ + 1 - 2k ) + \ (! - u ) - 2^) 
fe=i 

One can separate the real and imaginary parts of this expression. This helps in de- 
termining whether we require additional constraints on rapidities \j to ensure reality of 
the total energy. 

Numerical study showed that there exist several spurious solutions and only a subset 
of them, corresponding to a specific choice of counting numbers Ia \ is admissible. Non- 
string solutions exist, however they are not as numerous. From the numerical study we 
make the assumption that the elementary excitations over the FM (a sea of (2s)-strings 
) are (a) a pair of (1+) strings , and (b) (1+) and (1 — ) strings. 

The FM ground state is a filled band of (2s) strings. The density of ground state 
energy is 

N-l 

60 = mToo M = J d X/Wx)WA) - 2 E cot [jfj (78) 

The observed correlation between the integers suggest that the rapidities correspond- 
ing to (2s) and (a) ought to be connected allowing cancellation of the imaginary part of 
the total energy. 

Im [e (2s) ] = Im [e (a) ] (79) 
It can be shown that, Im [e( 2s )] = Im [ e (a)] an d Re [e(2«)] = —Re [e( n )] • 

& (80) 

Similar argument holds for (b)-type excitations, where Im [e(2 S )] = Im [ e (b)] an d Re [e(2s)] 
-Re [e (6) ]. 

One should note that the dressed energy and bare energy are equal since the function 
coupling the ground state density to excited states 6^. ^ is zero for j = (2s). Thus we 
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arrive at 

M («0 M (b) 

E = E q +Y: 2e (o) (x?) + £ 2e W (x?>) (82) 

0=1 0=1 

where E is the ground state energy. 

We now turn to the calculation of momentum. The momentum associated with a 
string of length j and parity v is found to be 

1~(D 



P(i+)(X) = -2 6 (1+)(X) 



1 



PO)(x) = -^eg( X ) /or j^(l+) (83) 
whence 

P(o)(x) = P(6)(x) = 2 arctan (tanh + tt (84) 

These expressions are similar to the ones for non-string excitations in the Anti- 
Ferromagnetic case for odd spin FZM. We get the dispersion relations for elementary 
excitations over the Ferromagnetic ground state as, 

C(0)(P) = 4sin (f) 

e (6) (p)=4sin(|) (85) 



VIII. Free Energy in the Anti-Ferromagnetic case for 
N even 

From the numerical study for finite lattices it was apparent that the Anti-Ferromagnetic 
(AFM) ground-state corresponds to a filled band of real roots for both T q and Tjj q . In 
other words the AFM ground state is a filled sea of (1,+) strings for both families of 
transfer matrices. 



The AFM ground state corresponds to real roots for both families of transfer matrices. 
Hence we consider the natural logarithm of both sides of BAE-s. From BAE(l) we get, 

vf i In ( ^MXj-Xk-^ =2M . n ( ^l±m + 2nIj m 
with standard definition of pi(x) and P2(x)j 

M 

2M-K 



Mx«) = !ef»'( X( »)-i-^ef>'( Xo -xrf 



0=1 



1 / 1 r-roo f 

From BAE(2) we get, 



(87) 



^ Vsinh(x,-Xfc + 2^7)y 1 
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In the continuum limit M — > oo we get 



P2(x) = (xa ~ M)pi(M) (89) 

The above pair of BAE is solved as before by the Fourier transform method. 

/ g l(fc) = IeW'(fc)-^-ef ) '(fc)p 2 (fc) (90) 

Mk) = (fe)pi(fc) (91) 



From the above two equations, 



sinh(^) . S jnh(^ - f|) 
sinh(^) • sinh(7rfc + 



Following a similar procedure as shown in detail in the Ferromagnetic case, the free 
energy for the Antifcrromagnctic case is obtained as, 

f = 7; I ~, — w ,\ cosh fc(?r - 4u - — ) - cosh fc(7T - — (93 
2J_ 00 k sinh(Trfc) L v N J v N'i y ' 

Substituting for p\{2k) 

. 1 f +00 dk sinh(Trfc-f) 7T. 

f = - -n- — ^-t— • cosh k(ir - 4m - — ) - 

J 2i_ oc fc sinh(^) ■ sinh(27rfc + ^) V AT 

sinh(7rfc-^) , 7T . 

/ -; To t s , 9 ^ ' cosh H^-lTr) (94) 

./_«, fc sinh(^) • sinh(27rfc + ^) V N J y ' 

IX. Excitation on AFM ground state 

In the case of AFM, the ground state is a filled band of (1+) . The results of previous 
section suggest that excitations appear as (2s) strings when a hole is created in the (1+) 
sea. Consider the Z(i+) (x) function. In the general case it should look like 

1 M <fc > 

%+)(*) = w - E £ e (£,*)(x - ( 95 ) 

fc /3=1 

The density of (1+) vacancies is given by 

<WX) = -^(1 +) (X) (96) 
The vacancy density <J(i+)(x) an d the density of (1+) particles P(i+)(x) 1S related by 

M< 1+) 

*(i+)(x) = P(i+)(x) + E 5 (x - X { p +)h ) (97) 

0=1 

where X^ 1+ ^ are the position of the holes. 
Thus 

M< fc > 

-Wx) = ^<)(x)-^ E E (L)(X-X^)- 

fc/(l+) 0=1 

M< 1+) 

^ / <,!+)(* - ^ + ^ E <,i + )(* - X? +)/l ) (98) 

J /3=1 



22 



The above equation can be interpreted as a collection of terms contributing to (1+) 
ground state, (1+) holes and excited particles over the ground state sea of (1+) string. 

C(i+) = °"(i+) + (7 (i+) + J2j CT (i+)' wnere CT (i+) i s the same as P(i+) of the last section. 

We obtain the energy of a state designated by a given set of strings {k}, having 
strings of type (fc) with string centers at Xp*\ 

e = E E^ fe) ) 

k 0=1 
strings 

M h +) M (fc) 

= / d X a (1+) (x)e il+) (x) - E e U+) {^ +)h ) + E E e (X?) 

J 0=1 + ) /3=1 

From the numerical study we can make the assumption that the elementary excita- 
tions over the AFM ground state, which is a sea of (l+)-strings, is given by a set of (2s) 
strings. 

The density of ground state energy for the AFM case is given by, 

N-l 

eo = J™o§ = I d xP(i +) (xhi +) (A) - 2 E cot (jf) (") 

Total energy is real since the imaginary part of the energy contribution from holes 
cancel the imaginary part of the energy of excitation, i.e., Ira [e(2s)] = I m [ e (i+)] ■ The 
real parts are given by, 

Re[e (2s) (x)]=Re[e {1+)h ( X )} (100) 
In this case also, the dressed energy equals the bare energy. As before, we will denote by 
£( 2s ) (xjj S ^ its real part. The total energy is given by, 

M (2s) 

E = E + E 2£ (2s) (X^) (101) 
0=1 

where £0 is the ground state energy. 

We now turn to the calculation of momentum. The momentum associated with a string 
of length j and parity v is found to be, 

f(i+)W = 4 e (!l)M ( 102 ) 

PU)(X) = -W ( )](X), forj ^(1+) (103) 



whence 

P(2s)(x) = 2arctan(tanh (|)) (104) 

If one keeps in mind the fact that the correlation of X s demand an additional tt for the 
creation of (1+) hole, we can derive the dispersion relation 

e (2s) (f>)=4sin(|) (105) 
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